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Abstract 

The forward limit of the chiral-odd generalized parton distributions (GPDs) and their lower 
moments are investigated within the framework of the chiral quark soliton model (CQSM), with 
particular emphasis upon the transversity decomposition of nucleon angular momentum proposed 
by Burkardt. A strong correlation between quark spin and orbital angular momentum inside the 
nucleon is manifest itself in the derived second moment sum rule within the CQSM, thereby pro- 
viding with an additional support to the qualitative connection between chiral-odd GPDs and the 
Boer-Mulders effects. We further confirm isoscalar dominance of the corresponding first moment 
sum rule, which indicates that the Boer-Mulders functions for the u- and d-quarks have roughly 
equal magnitude with the same sign. Also made are some comments on the recent empirical extrac- 
tion of the tensor charges of the nucleon by Anselmino et al. We demonstrate that a comparison 
of their result with any theoretical predictions must be done with great care, in consideration of 
fairly strong scale dependence of tensor charges, especially at lower renormalization scale. 
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I. INTRODUCTION 



The concept of generalized parton distributions (GPDs) has recently attracted consid- 
erable interest Naturally, these new quantities contain richer information on the 
internal quark-gluon structure of the nucleon, well beyond what can be learned from the 
usual parton distribution functions (PDFs). A complete set of quark GPDs at the leading 
twist 2 contains four helicity conserving distributions, usually denoted as H q , E q , H q , E q , 
and four helicity-fhp (chiral-odd) distributions, labeled as Hj,,E^,H^,E^ These 
GPDs are all functions of three kinematical variables, x, £, and t, where a; is a generalized 
Bjorken variable, t is the four-momentum-transfer square of the nucleon, while £ is the lon- 
gitudinal momentum transfer, usually called the skewdness parameter. The standard PFDs 
are naturally contained as a subset of these GPDs. That is, in the forward limit £, t — > of 
zero momentum transfer, H q (x,£,t) and H q (x,£,t) reduce to the unpolarized distribution 
function q(x), and the longitudinally polarized distribution functions Aq(x), respectively. 
On the other hand, Hj>(x,£,t) reduces to the so-called transversity distribution function 
A T q(x). 

Experimental studies so far have mostly been concentrated on the helicity-conserving 
(chiral-even) GPDs, especially on H q (x,^,t) and E q (x,£,t) Q| -12], because they are very 



interesting quantities for clarifying the role of quark orbital angular momentum in the nu- 
cleon spin problem and also because they are easier to access experimentally 
as compared with the helicity-flip (chiral-odd) GPDs. Although there exist some proposals 
to access to the chiral-odd GPDs in diffractive double meson production 161 ]. 171 ]. we now 
have almost no empirical information on them. (An exception is the forward limit of GPD 
H q (x,£,t), i.e. the transversity Aj<q(x). The first empirical extraction of the transversity 
distribution has recently been done by Anselmino et al. based on the combined global 
analysis of the measured azimuthal asymmetries in semi-inclusive scatterings and those in 
e + e~ — ► hih 2 X processes fl^.[l9].) 

Although a direct experimental access to the chiral-odd GPDs is not very easy at the 
present moment, it was shown by Burkardt that they are not only interesting from a theoret- 
ical viewpoint but also they have important influence on some physical observables [20!]. [21]. 
First, the transversity decomposition of quark angular momentum in the nucleon intro- 
duced by him indicates that a strong correlation between quark spin and angular momen- 
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turn is hidden in the 2nd moment of chiral-odd GPDs, more specifically in the combination 
iif + 2 iff. + Ej,. He also suggested that a strong correlation would exist between the 1st 
moment of 2 H q + E q and the Boer-Mulders functions h\ q describing the asymmetry of the 
transverse momentum of quarks perpendicular to the quark spin in an unpolarized target 
22]. (This is a variant of the analogous relation between the Sivers function y and the 



mm-, 



anomalous magnetic moment of a quark with the flavor q also proposed by him 

Turning to the status of theoretical studies of chiral-odd GPDs, most works so far have 
been restricted to the studies of the transversity Aq(x), which is the forward limit of the 
GPD Hj<(x,£,t). A lot of model calculations were reported on the transversity and its 1st 
moment, i.e. the tensor charge 26] - 37J. There also exist lattice QCD studies on the lower 
moments of i?f.(x, £, t). Its 1st moment, i.e. the tensor charge, was first investigated in 38], 
while the simulations were extended to include its 2nd moment as well in [39|. However, 
the lattice QCD studies on the moments of other GPDs, i.e. Ej,, iff, E\ , have not been 
reported yet. Concerning the full x, £, t dependence of the chiral-odd GPDs, there have 



been only a 



'ew model calculations. The one is the investigation by Pasquini, Pincetti, and 



Boffi 401 ]- 42j within the framework of the light- 



similar investigation by Dahiya and Mukherjee 



ront constituent quark model (see also a 



43]), and the other is the calculation by 



Scopetta based on a simple version of the MIT bag model. Probably, most extensive is the 



investigation by Pasquini et al. 



40J . They gave predictions not only for the lower moments 



of GPDs but also for the full x, £, t-dependence of those GPDs. Note however that their 
calculations were made possible in price of one crude approximation. That is, in their model 
calculations, only the lowest-order Fock-space components of the light-front wave functions 
with three valence quarks are taken into account. The previous investigation of the chiral- 
even unpolarized GPD ff( /=0 )(x, £, t) based on the approximate treatment of the CQSM 



45] indicates that this approximation is not necessarily justified, and the inclusion of higher 



Fock- components may bring about richer x- and ^-dependence in the GPDs. 

In the present investigation, we try to investigate chiral-odd GPDs beyond the three 
valence quark approximation. Within the CQSM, which we shall use, the effects of higher 
Fock-space components can be included nonperturbatively as contributions of deformed 



Dirac-sea quarks in the hedgehog mean field 



52], [53], not through the perturbative Fock- 



space expansion. Unfortunately, technical hardness of this ambitious program does not allow 
us complete calculation of GPDs in full dependence of the three kinematical variables, x, £ 



3 



and t, at the present stage. In the present investigation, we therefore content ourselves in 



the calculation of the forward limit of a GPD, i.e. Gj<(x, 0, 0) = lim^_ 



•o,t->o 



2H^(x, £,,t) + Erp(x,£,t)], since, except for the transversity distribution H^(x, 0,0), this 
distribution is physically the most interesting quantity, which is thought to contain valuable 
information on the correlation between the quark spin and orbital angular momentum inside 
the nucleon as pointed out by Burkardt [^(J, {21]. (There already exist several investigations 
within the CQSM on the forward limit of chiral-even GPD E(x, £ , t) 
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471 as well as on 
5l|.) Concerning the 



the generalized form factors corresponding to its lower moments 
transversity Atq(x) = H^(x, 0, 0) contained in the combination to define G^(x, 0, 0), the first 
empirical information has recently obtained by Anselmino et al. through the semi-inclusive 



deep inelastic scatterings 



18] 1 ISfl- The extracted transversities and/or their 1st moments 



were compared with some model predictions. We shall discuss in the present paper, that 
such a comparison is potentially very dangerous if one does not pay the closest attention to 
the fairly strong scale dependence of the transversities. 

The paper is organized as follows. First, in sect Jill we shall derive theoretical formulas, 
which are necessary for evaluating the relevant GPDs within the framework of the CQSM. 
Next, in the first part of sect JUIj we show the results of numerical calculation for the isoscalar 
and isovector part of G^{x, 0, 0) as well as their 1st and 2nd moments. The second part of 
sect JIIII is devoted to the discussion on the delicacy, which is shown to arise in a comparison 
between the recent empirical determination of the tensor charges and the corresponding 
theoretical predictions. Some concluding remarks are then given in sect HVI 

II. CHIRAL-ODD GPDS IN THE CQSM 



The chiral-odd GPDs are defined as nonforward matrix elements of light-cone correlation 
of the tensor current as 



M 



2P+ 



u(p',X' 



H T (x,£,t)ia + i l5 + H T (x,£,t) 



Ml 



j f +j°<i3 A ry. j f +ja/3 p ^ 

+ E T (x,£,t) — + Er(x,£,t)- 



2M 
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M 
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u(p,\), (1) 
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where % = 1, 2 is a transverse index, while p (p') and A (A') are the momentum and the 
helicity of the initial (final) nucleon, respectively. We use here the light-cone coordinates 



v = [v 



± v 3 ) I y/2, and v± = (v 1 , v 2 ) for any four- vector v^. We also use the notation 



P = + P)> A = P'-P- 



(2) 



As is widely known, the GPDs depends on three kinematical variables, x, £, and t, where 
x is a generalized Bjorken variable, t = A 2 is the four-momentum-transfer squared of the 
nucleon, and £ = — A + / (2 P + ) denotes the longitudinal momentum transfer, usually called 
the skewdness parameter. 

For model calculation, it is convenient to work in the so-called Breit frame, in which 



P 



E A/2 , + A/2), 



P 



E 



A/ 2 ,- A /2 



(3) 



so that 



P = (J5 A/2J 0), A = (0,A). 



(4) 



We also assume large N c kinematics, in which the nucleon is heavy, ~ 0(N C ), and its 
center-of-mass motion is essentially nonrelativistic. Under these circumstances, A 1 = O(iVj) 
and A = C^iV^ 1 ), so that the hierarchy holds that M N > |A*| > |A°|. Note also that t = 
- A = 0(iV c °) and £ = - A 3 /(2Mjv) = 0(N C 7 1 ). Then, noting that A = (A ± ,-2M N £), 
we evaluate the right hand side of eq. (pQ), to obtain 



M 



H T + 



8M 2 N 



E T — - H T ) - £ E T 
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2 

1 

2M N 
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(5) 



Since we are interested in the forward limit (Aj_ — * 0, £ — > 0) in the present investigation, 
we can project out these four independent pieces as 
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M, (9) 



where is the azimuthal angle of the transverse vector Aj_, i.e. A_ 
For convenience, let us use below the shorthand notation : 



. cos 0, sin i 



G T {x,£,t) = H T (x,£,t) + 2H T (x,£,t) + E T (x,£,t), 
K T (x,£,t) = Er(x,Z,t) - ^H T (x,£,t). 



(10) 
(11) 



Now, we are ready to evaluate the amplitudes Ai explicitly in the CQSM. We first 
investigate the answer at the mean field level, i.e. we derive theoretical expressions for 
the 0(Q°) contribution to A4, with Q being the collective angular velocity of the rotating 
ledgehog mean field. Using the formalism developed in the previous studies 32|, 54j - 
6ll |. the isoscalar part of A4(x, 0, t) at the 0(Q°) is given as a sum over all the occupied 
eigen-states of the Dirac Hamiltonian H : 

dzQ 



M {I=0) (x,0,t) = M N j — / d 6 x e 

2 7T 



3^, „iA.±-X 



X 



N c £ e i*o(*M N -E n ) $t( :r )( 7l75 _, 72 )$ n ( a ,_ 2oe 3) 



;i2) 



ndocc 



with e3 = (0, 0, 1) being a unit vector in the z direction. Here, $ n (JC) are the eigenfunctions 
of the Dirac hamiltonian with the hedgehog mean field, i.e. 



H$ n (x) = E n <f> n (x) 



(13) 



with 



Noting that 



we easily find that 



i4 7 - 0) (x,0,0) = 



,%A ± -x 



1 + i Aj_ ■ x - - (Aj_ • x) 2 + 



G^- 0) (x,0,0) = M N N C J] (n\x 2 ( lll5 -i l2 )S(xM N -E n -p 3 )\n), 



2M, 



N 



(14) 
(15) 

(16) 
(17) 



ndocc 
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2M 



' ^ /=0) (x,0,0) = 0, (18) 



N 

^K { T I=0 \x,0,0) = 0, (19) 



4M N 



which shows that only G^~°\x, 0, 0) survives among the four GPDs, at the lowest order in 
Q, or in 1/N C expansion. 

Next, we turn to the isovector part. The isovector part of Ai(x, 0, t) is given as 

M {I=1 \x } 0,t) = M N [ — [ d 3 xe iA ^ x 

J 2 7T J 

N c J2 e**^""-^ & n {x) A* r 3 A( 7l75 - i l2 )® n {x - z e 3 ), (20) 



x 

ngocc 



where A is the rotation matrix belonging to flavor SU(2). Here we use the identity 

A*r 3 A = D 3k (A)r k , (21) 

where D 3k (A) is a Wigner's rotation matrix, which should eventually be sandwiched be- 
tween the rotational wave functions *&t 3 ~j 3 ~ \A\ representing the spin-isospin states of the 
nucleon. Using the expansion (fl5l) . together with the replacement 

D 3k {A) — -+ - - (r 3 ) T , Ta {a k ) s , S3 , (22) 

which should be interpreted as an abbreviation of the identity 

/ 4t4 =1/2) *H D ^ A ) *S =1/2) H VA = ~\ ^t>t 3 Ms>s 3 . (23) 

in the projection formulas ([7]) - (J9j), we are led to the following expressions for the O(Q ) 
contributions to the isovector parts of four GPDs as 

i4 /=1) (x,0,0) = M N (-^) (n\r 1 ( lll5 -t l2 )5(xM N -E n -p 3 )\n), (24) 



n£occ 



2M 



' G? _1) (x,0,0) = 0, (25) 



N 



1 4 /_1) (^0,0) = M N £ ( n I ' x i T 3 (7i 75 - * 72) S(xM N -E n - p 3 ) I n), (26) 



2M N \ 3 / ngocc 

1 „rr=n, . / N, 



x (n I - [2 (x± ■ t)xi - x]_ti] (7175 - i^ 2 ) S(xM N - E n - p 3 ) \n). (27) 

Note that, just opposite to the isoscalar case, only G ( ^~ 1 \x, 0, 0) vanishes, while other three 
GPDs are generally nonzero. 



The O^ 1 ) contributions, or equivalently, the next-to-leading contributions in 1/N C ex- 
pansion, can similarly be evaluated, although the manipulation is much more complicated. 
Skipping the detailed derivation, we first write down the answers for the isoscalar parts : 

4 J=0) (x,0,0) = — M N Y, (m\n\n) 

m£all,n£occ 

x (n | ( 71 75 - i 72 ) [^-L- 6(x M N -E n - P3 ) | m) , (28) 

' G^ =0) (x,0,0) = 0, (29) 



2M N 

' 4 /=0) (x,0,0) = -M N ^- J2 {m\T 3 \n 



2M N 1 v ' ' ' 21 c 

x (n\ix 1 (71 75-^72) (^4^-^^) <H* M " - K - Pa) | m), (30) 
^JfCf^foO.O) = — M N J2 (m\r c \n) 

N m£all,n£occ 

x ( n I g t ( x i ~~ ^a) &=,i + 2 xi x 2 5 C , 2 ] (71 75 - i 72) 

* (^-2^^)^ M "- B "-^' m >' (31) 
On the other hand, the O^Q 1 ) contributions to the isovector part are given as 



^- ij (x,0,0) = te lac M N -j £ 



fi /" ^ E — E 



x (m I r c I n) (ra I 7^(7175 - 173) ^(icMjv - - p 3 ) \ m), (32) 



1 4 /=1) (x,0,0) = ~M N ^- £ ("»kc|n> 



2Mjv J 6/ 



rndall ,n£occ 



x (n|r c a; 2 ( 7 i 7 5-i 72 ) L__ - _L_ A^J 5(x Mjv - E n - p 3 ) | m), (33) 

-i-4^^,0,0) = .. 3ac M w § E Tr^r 

^ iw iV u J m£nocc,n£occ ^rn 

x (m| r c \ n) (n \ r a ix 1 ( 7l 75 -in/2) S(x M N - E n - p 3 ) \ m), (34) 
— ^-^ /=1) (x,0,0) = ie bac M N -j E 7-1 1 7-1 M^c lw) 

x ( n I r a ^ [ ki - ^2) <W + 2x 2 a; 2 (5 fei2 ] ( 7 i7 5 - 272)5(0;^^^ - E n -p 3 ) |m) . (35) 

To sum up, we can summarize the novel f2 dependence of the eight GPDs as 

H^ I=o \x,0,0) = + OiQ 1 ), (36) 
G^ =0) (a;,0,0) = O(Q ) + 0, (37) 
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£^- 0) (x,0,0) = + O^ 1 ), (38) 

K^ I=0 \x,0,0) = + OiQ 1 ), (39) 

and 

^ /=1) (a;,0,0) = O(fi ) + O^ 1 ), (40) 

G^ /=1) (x,0,0) = + OiQ 1 ), (41) 

i^ /=1) (x,0,0) = 0{Q°) + OiQ 1 ), (42) 

K { T I=1 \x,0,0) = 0(n°) + O^ 1 ), (43) 

where terms which do not contribute are denoted as 0. Since Q is an 0(1/N C ) quantity, 
this especially means that Gt(x, 0, 0) is a quantity with isoscalar dominance, although the 
isovector component also exists as an 1/N C correction. (See the discussion in sect. III.) One 
may also notice that the contributions to the isovector GPDs, H? 1 ', E^ 1 ', and , 
survive at the mean-field level, or at the O(Q ) level, while these GPDs receive O^ 1 ) 
contributions as well. The appearance of the antisymmetric e-tensor, as observed in Eqs. (l40l) . 
( 142]) and ( l4"3"l) , is a characteristic feature of this novel 1/N C correction. This unique 1/N C 
correction is known to play an important role for resolving the notorious underestimation 
problem of some isovector observables of the nucleon like the isovector axial charge and the 



62], leal, [64 1 



isovector magnetic moment, inherent in the hedgehog-type soliton model 

So far, we have derived the theoretical expressions for the forward limits of four chiral- 
odd GPDs, H T (x,£,t), E T (x,£,t), H T (x,£,t), and E T (x,£,t) in the CQSM. Since the for- 
ward limit of Hx(x,^,t), which is known to reduce to the familiar transversity distribution 



ATq(x), has already been investigated within the CQSM [32] -[36(|, we need to evaluate 
the remaining three independent GPDs, E T (x, 0,0), H T (x,0,0), and E T (x,0, 0), or equiv- 
alently, Gt(x, 0,0), K T (x,0,0), and E T (x,0,0). Unfortunately, we find that the numerical 
calculation of Kt(x, 0,0) is quite involved. In the present study, we therefore concentrate 
on Gt(x, 0,0), which is a special combination of three GPDs, -£^(^,0,0), Ht(x, 0,0), and 
Et{x, 0, 0), as given by ( FlOl) . From a physical viewpoint, this is the most interesting quan- 
tity, which appears in Burkardt's transversity decomposition of quark angular momentum 

3,0- 

According to Burkardt, the transverse decomposition of angular momentum in the nu- 
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cleon is given in the form 

<•/?> = (-W + (■£,-*>> (44) 
where the 1st and the 2nd terms in the right hand side respectively stand for the angular 
momentum carried by quarks with transverse polarization in the + x and — x directions in 
an unpolarized nucleon at rest. On the other hand, the difference of the above two quantities 
gives the transverse asymmetry 

(S x J x q ) = K+*> - (45) 

which can be interpreted as representing a correlation between quark spin and orbital angular 
momentum in an unpolarized nucleon. Burkardt has derived the identities, which relate the 
above two quantities to the 1st and the 2nd moment of GPDs as 

{J x q ) = y J\ [H(x,0,0) + E(x,0,0)} dx (46) 

{5J x q ) = \ J X x [H T (x,0,0) + 2H T (x,0,0) + E T (x,0,0) 
1 r 1 



dx 

, xG T (x,0,0)dx. (47) 
2 J — i 

The quantities appearing in the 1st sum rule are the forward limit of the familiar (chiral- 
even) unpolarized GPDs, H(x,£,t) and E(x,£,t). This identity is essentially Ji's nucleon 
spin sum rule and nothing new [3^]. What is new is the second sum rule. It relates the above- 
mentioned transverse asymmetry to the 2nd moment of the chiral-odd GPD Gt, which we 
recall is a particular combination of Ht, Ht, and Et- It is interesting to see the explicit 
form for the 2nd moment of G^~°\x, 0, 0) in the CQSM. From (JT7I) . we readily find that 



J i xG ( T =0 \x,0,0)dx 

N c \E n {n\{-i) 1 -x\n) + ]- (n | T ° S • L \ n) } 

n£occ •* 



2 
3 



2 ^ / -icr ■ x \ 1 cr ■ L 

N c E K(n\ \ in) + -(n| | n) , (48) 



J neocc \icr-x / z \ — cr ■ L 

with S = 7° 775 being the relativistic spin operator of quark field. As argued by Burkardt, 
the transverse asymmetry signals the correlation between quark spin and orbital angular 
momentum in an unpolarized target. One can clearly see that such correlation manifests 
itself in the 2nd term of the above sum rule (jUJ), since it reduces to the nucleon matrix 

10 



element (at the mean-field level) of the operator 7 X • L, which is certainly the scalar 
product of the relativistic spin and orbital angular momentum of quarks aside from an extra 
factor 7 . Unfortunately, the 1st term of eq. (l48p is a highly model dependent expression, as 
convinced from the appearance of the single-particle energy E n of the Dirac Hamiltonian H 
with the hedgehog mean field. This makes a simple physical interpretation of the 1st term 
not so easy. 

Before ending this section, we want to make a brief comment on the GPD ET{x,^,t). 
Although this GPD is generally nonzero (see Eqs. (126|) . (1341) . and (130|) ). its 1st moment is 
known to vanish by time reversal invariance Q. As a consistency check of our theoretical 
framework, we shall explicitly prove in Appendix that the 1st moment of Et(x, 0, 0) in fact 
vanishes identically. 



III. NUMERICAL RESULTS AND DISCUSSIONS 



A. chiral-odd GPDs and transversity decomposition of angular momentum 

Within the framework of the CQSM, the expression of any nucleon observable is divided 
into two parts, i.e. the contribution of what-we-call the valence quark level (it is the lowest 
energy eigenstate of a Dirac equation with the hedgehog mean field, which emerges from 
the positive energy continuum) and that of the deformed Dirac sea quarks. Since the latter 
contains ultraviolet divergences, it must be regularized. Here, we use the Pauli-Villars 



regularization scheme with single subtraction, for simplicity [54|-[32|]. The Pauli-Villars 
regulator mass Mpy is not an adjustable parameter of the model. It is uniquely determined 
from a model consistency, once the dynamical quark mass M, the only one parameter of 
the CQSM, is fixed to be M = 375 MeV from the phenomenology of the nucleon low energy 
observables. 

We first show in FigJT] the CQSM predictions for the forward limit of the isoscalar 
GPD G^ =0) (x,£,t), i.e. G^ =0) (x, 0, 0) = lim^ ,t-o [H^ =0 \x,^,t) + 2 H^ I=0) {x,^t) + 
E^~°\x,^,t)}. Here, the dashed and dash-dotted curves respectively stand for the contri- 
bution of N c ( — 3 ) valence quarks and that of deformed Dirac-sea quarks, while their sum 
is shown by the solid curve. The distribution function in the negative x region should be 
interpreted as antiquark distribution except for an extra minus sign related to the charge- 
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FIG. 1: The prediction of the CQSM for the forward limit of the isoscalar GPD Gj< (a;, 0,0) = 
lim£_>o,t->o 

[H^~°\x,C,t) +2H^-°\x,C,t) + E$-°\x,£,t)]. The dashed and dotted curves re- 
spectively stand for the contributions oi N c ( = 3) valence quarks and of deformed Dirac-sea quarks, 
while their sum is shown by the solid curve. 

conjugation property of this distribution. One clearly sees a strong chiral enhancement of the 
deformed Dirac-sea contribution in the small x region. We recall the fact that a similar chiral 
enhancement of the Dirac-sea contribution is also observed in the CQSM prediction for more 
familiar unpolarized parton distribution function of isoscalar type, and that it plays a crucial 
role for ensuring the positivity condition of the antiquark distribution u(x) + d(x) [541]. [55]. 
Naturally, such chiral enhancement of the antiquark distribution cannot be reproduced by 
such a model as light-cone constituent quark model with N c ( = 3) quark approximation. 
Next, shown in Figj2] is the CQSM prediction for the forward limit of the isovector 



GPD G^~ 1 \x,^,t), i.e. G? _1) (x, 0, 0) = lim e _ ,t-*o [H$- 1} (x,£,t) + 2H^- 1, (x,^t) + 
Ej. (x } £,£)]. Here, the meaning of the curves is the same as in the previous figure. Also for 
the isovector distribution, one observes a strong chiral enhancement of the deformed Dirac- 
sea contribution in the small x region. The x dependence of the Dirac-sea contribution 
for this isovector distribution turns out to be totally different from the isoscalar distribu- 



r('=i)/ 



r('=i)/ 
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tion, however. The deformed Dirac-sea contribution for the isovector distribution is nearly 
symmetric with respect to the variable change x — > — x, in sharp contrast to the isoscalar 
distribution, which is approximately antisymmetric. This behavior i s ag ain resembling more 



familiar unpolarized parton distribution function of isovector type [58|,[59(, 60J. We recall 
the fact that this chiral enhancement of the isovector unpolarized distribution is just what 
is required by the celebrated NMC measurement, which established the dominance of d sea 



over the u sea inside the proton 



651 ] . Unfortunately, we do not have any simple explanation 



about why such a similarity exists between the small x behaviors of the unpolarized parton 
distribution functions q(x) and the forward limit of the GPD Gt(x, £, t). 



G { t =1] (x,0,0) 




-1 1 — ■ — ^ — ■ — ■ — ' — ^ — ■ — ■ — — 1 — — — ^ — — — ^ — ■ — ■ — 1 
-1.0 -0.5 0.0 0.5 1.0 

x 

FIG. 2: The prediction of the CQSM for the forward limit of the isovector GPD G^ =1) (x,0,0) = 
lim^o.t^o [H^ (x,£,t) + 2H£ 1 \x,£,t) + E^~ 1 \x, £, t) ]. The meaning of the curves is the 
same as in FigfTj 

As discussed before, the above two distribution, i.e, G^~°\x, 0,0) and (x, 0, 0), 

are very interesting quantities from a physical viewpoint, since they give the transversity 
decomposition of the angular momentum inside the nucleon. So far, there exist only a few 
theoretical investigations on the chiral-odd GPDs and their moments. In table [J we com- 
pare the CQSM predictions for the transverse asymmetries with those of the two versions 
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of the light-front constituent quark model by Pasquini et al., that is the harmonic oscilla- 



tor (HO) model and the hyper central model 40j. Note that their models are essentially 
the three quark model with relativistic kinematics. One finds that the predictions of the 
CQSM just lie between the HO model and the hyper central model. In all the three mod- 
els, the isoscalar transverse asymmetry is seen to be larger than the isovector one, but the 
isoscalar-to-isovector ratio is largest in the CQSM. This observation is qualitatively consis- 



tent with the large N c prediction given in 66]. However, literally taking the N c — > oo limit, 
the ratio (S X J^ ) / (8 X J% ) would become infinite. Our present analysis here shows that 
GT(x,£,t) is an isoscalar-dominant quantity but the isovector component, which arises as a 
1/N C correction, is also important. 



transverse asymmetry 


HO 


Hypercentral 


CQSM 


(8*JX) 


0.68 


0.39 


0.49 




0.28 


0.10 


0.22 


2(S x J% +d ) 


1.92 


0.98 


1.41 


2{8*JT d ) 


0.80 


0.58 


0.54 


(5*J^ d )/(5*JZ- d ) 


2.40 


1.69 


2.61 



TABLE I: Some theoretical predictions for the transverse asymmetry. Here, the second and the 
third columns stand for the two versions of the light-front constituent quark model of Pasquini et 
al. 40j], i.e. the harmonic oscillator model (HO) and the hyper central model (HYP), while the 



predictions of the CQSM are shown in the fourth column. 

Also interesting is the 1st moment sum rule for Gt, which can be divided into two 
pieces, i.e. the 1st moment of the transversity H^(x, 0, 0) and that of the distribution 
2if|(x,0,0) + E^(x,0,0) as 

j G q T (x,0,0)dx = J 1 H^(x,0,0)dx + J 1 [2^(x,0,0) + E%(x,0,0)]dx. (49) 

Here, the 1st term of the r.h.s. of the above equation, i.e. the 1st moment of the transversity, 
gives the tensor charge 

A T q = f 1 H^(x,0,0)dx. (50) 
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On the other hand, the 2nd term, i.e. the 1st moment of 2 H T (x, 0, 0) + E T (x, 0, 0) defined 
by 

4 = J [2 H T (x, 0,0) + E T (x,0,0)]dx, (51) 

was given an interpretation as a quantity governing the transverse spin-flavor dipole moment 
in an unpolarized target by Burkardt. In fact, he showed that k t gives us an information on 
how far and in which direction the average position of quarks with spin in the x direction 
for an unpolarized target relative to the center of momentum. The decomposition (149p 
corresponds to a similar decomposition of the 1st moment sum rule for the unpolarized GPD, 
Eli(x,^,t) = H q (x,£,t) + E q (x, £, t), which gives the total magnetic moment consisting of 
the quark number N q and the anomalous magnetic moment n q as 

J i E q M (x, 0, 0) dx = J^H q (x,0,0)dx + J ^ E q (x, 0, 0) dx = N q + n q . (52) 

In tableHIl we again compare the CQSM predictions for the quantities k t (here we tentatively 
call it the "anomalous tensor moment") with the corresponding predictions of Pasquini et 
al. Here, the prediction for the isoscalar part is closer to that of the HO model, while the 
prediction for the isovector part is closer to that of the hyper central model. 



1st moment of 2 Ht + Et 


HO 


Hypercentral 


CQSM 


FXirp 


3.60 


1.98 


3.47 


t\jrp 


2.36 


1.17 


2.60 


K u+d 

IX rp 


5.96 


3.15 


6.07 


K u-d 


1.24 


0.81 


0.88 


u+d / u-d 

r\jrp J r\,rp 


4.81 


3.89 


6.90 



TABLE II: The theoretical predictions for the 1st moment of 2Ht + Et- 



According to Burkardt's conjecture, one would expect an intimate connection between 
the time-reversal-odd (T-odd) transverse momentum-dependent distributions and the GPDs. 
They are the approximate proportionality relation between Siver's function and the anoma- 
lous magnetic moment with opposite sign 24|, 251 ]. 



AC 9 , 



(53) 
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and also the proportionality relation between Boer-Mulders' function and the anomalous 



tensor moment 



mm, 

hi q (x,k 2 ± ) 4. (54) 



If his conjecture is combined with with some typical model predictions for the anomalous 
tensor moments, one would get the following approximate relations : 

h^ d ~ -hi u : MIT bag model, 
h± d ~ h± u : Large N c prediction, 
h± d ~ -h± v : CQSM, 

thereby dictating that the Boer-Mulders functions for the it- and ci-quarks would have the 
same sign, although the predictions on the relative magnitudes are a little variant. This 
should be contrasted with the fact that Sivers functions for the u- and <i-quarks appears to 
have opposite sign as 

fi d ~ ~fi u , (55) 
in conformity with empirically known relation 

n d k u . (56) 

From our viewpoint, the origin of this qualitative difference is very simple. It comes from the 
fact that the anomalous magnetic moment is a quantity with isovector dominance, whereas 
the quantities k% is of isoscalar dominance, as expected from the N c counting rule indicated 
in Eqs.flMD-fllS])- 



B. tensor charges : current empirical information versus theoretical predictions 

Some years ago, the first empirical extraction of the transversity distributions has been 
made by Anselmino et al. based on the combined global analysis of the measured azimuthal 
asymmetries in semi-inclusive deep inelastic scatterings (SID IS) and those in e + e~ — > h\h 2 X 



processes. More recently, they have further refined their g' 
from HERMES, COMPASS, and BELLE Collaborations 



obal analysis by using new data 



3- 



The 1st x- moments of the 



transversity distributions - related to the tensor charge - have been extracted to be 

A T u = 0.59 + °" 14 , A T d = -0.20 + °'° 5 , (57) 
- 0.13 - 0.07 
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at the renormalization scale Q 2 = 0.8 GeV 2 . They concluded that their new transversity 

distributions are close to some model predictions, especially the predictions by a covariant 

quark-diquark model by Cloet at al. [37]. This agreement is related to the fact that the 

predictions by Cloet et al. give the smallest magnitudes of tensor charges among many 

theoretical predictions including those of the lattice QCD. As we shall discuss below, this 

statement appears very misleading, however. A delicate point is that the tensor charges 

are strongly scale-dependent quantities especially at low renormalization scale. In fact, the 

I I 

bare predictions of the covariant quark-diquark model given in [37J for the tensor charges 
are nothing small. They are 

A T u = 1.04, A T d = -0.24, (58) 

or in the isospin language, 

A T g (/=1) = 1.28, A T g (/=0) = 0.80. (59) 

Cloet et al. regard the transversity distributions, which give the above 1st moments, as 
initial distributions given at the scale Q 2 = 0.16 GeV 2 , and take account of their scale 
dependencies by using the next-to-leading (NLO) evolution equation. This procedure gives 
the tensor charges at the scale Q 2 = 0.4 GeV 2 : 

A T u = 0.69, A T d = -0.16, (60) 

or equivalent ly 

A T q {I=1) = 0.85, A T q {I=0) = 0.53, (61) 

which are much smaller than the bare predictions of the model despite pretty small scale 
difference. 

As naturally anticipated, to start the NLO evolution at such low energy scale as Q 2 = 
0.16 GeV 2 is very dangerous. To convince it more concretely, we first show in Figj3]the QCD 
running coupling constant as(Q 2 ) at the NLO as a function of Q 2 . Here, we have used the 
standard NLO formula 



a s (Q ) 



47T 



P HQ 2 /A 2 ) 
where 



Pi lnln(Q 2 /A 2 ^ 



1 - 



P 2 ln(Q 2 /A 2 ) 



(62) 



2 38 
Po = 11 - -n fl fix = 102 - —n f , (63) 
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ols{Q 2 ) at next-to-leading order 



2.0 
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0.0 

12 3 4 

Q 2 [GeV 2 ] 

FIG. 3: The QCD running coupling constant as(Q 2 ) at the NLO in dependence of Q 2 , obtained 
with the effective flavor number nj = 3 and the QCD scale parameter A-^j^ = 0.248 GeV. The 
energy scale Q 2 = 0.16 GeV 2 and Q 2 = 0.40 GeV 2 are marked by open squares as a guide. 




together with the effective flavor number rif = 3 and the QCD scale parameter A = A.jjg = 
0.248 GeV, taken from the NLO analysis by Gliick, Reya, and Vogt [67]]. One sees that, at 
Q 2 = 0.16 GeV 2 , as(Q 2 ) is about 1.5, which immediately throw doubt on the use of the 
perturbative QCD evolution equation. 

The statement can be made more explicit by investigating the NLO evolution of the 
tensor charges themselves. The anomalous dimensions at the NLO, which control the scale 
dependencies of the moments of the transversities are given in 68], 69], 70j. We are interested 
here in the NLO evolution of the 1st moment, i.e. the tensor charges. (Note that, since the 
transversities do not couple to the gluon distributions, the evolution of the tensor charges 
is flavor independent. For more detail, see the discussion later.) The solution of the NLO 
evolution equation for the tensor charge A T q(Q 2 ) is given as 



A T q(Q 2 
A T q(fi 2 



' as(Q 2 
s a s (fi 2 



7 (0) /2/?o 



A) + /WQ 2 )/4tt 
Po + /?i«5(/i 2 )/4vr 



i^W/ft - 7 C°)//3 ) 



(64) 
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with 



7 (1) m 

7 (0) /2/3o 



724 



104 



n 



9 27 
4/(33-2n/) 



/ 2 102 



38 



■ n 



f ' 



(65) 
(66) 



To NLO accuracy, the above solution are sometimes expanded as 



A T q(Q 2 



« S (Q 2 



Y (0) /2/3o 



"s(/i 2 ) / 

'*smV /27 



1 _ J_& 

4vr/5 



7 



(i) 



7 



(o) 



2/5i 2/5 



«5(/i 2 ) -«5(Q 2 )] 



1 - 



337 
486 n 



K(/i 2 ) - « 5 (g 2 ) 



(67) 



Here, we have set rif = 3, which reproduces the form used in 37J. For large enough Q 2 , 
where the QCD running coupling constant is much smaller than unity, both expressions 
should approximately be equivalent. However, we have already pointed out that, at the 
scale of Q 2 = 0.16 GeV 2 , as is even larger than 1.5. 



A T q(Q 2 )/A T q(Q 2 im ) at NLO with Q 2 ni = 0.16 GeV 2 
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FIG. 4: The scale dependence of the tensor charge, where the evolution is started at fi 2 = Q 2 ni = 
0.16 GeV 2 . The solid and dashed curves respectively correspond to the results obtained with the 
exact (Eq. (|64p ) and approximate (Eq. (|67p solution of the NLO evolution equation. 

Shown in FigJH are the Q 2 -dependence of tensor charge, in which the evolution is started 
at /i 2 = Q 2 ni = 0.16 GeV 2 . The solid and dashed curves respectively correspond to the 
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answers obtained by using the exact (Eq. (l64j) ) and approximate (Eq. fl67|) ) solutions of the 
NLO evolution equation. One clearly observes a drastic difference between the two choices. 

A T q(Q 2 )/A T q(Ql t ) at NLO with Q 2 m = 0.34 GeV 2 

1.2 
1.1 
1.0 
0.9 
0.8 
0.7 
0.6 
0.5 

12 3 4 

Q 2 [GeV 2 ] 

FIG. 5: The scale dependence of the tensor charge, where the evolution is started at fj, 2 = Q 2 ni = 
0.34 GeV 2 . The solid and dashed curves respectively correspond to the results obtained with the 
exact (Eq. ([64"|) ) and approximate (Eq. ([6T)) ) solution of the NLO evolution equation. 

On the other hand, shown in FigfS] are the Q 2 dependence of the same quantity, where 
the evolution is started at /i 2 = Q 2 ni = 0.34 GeV 2 , which corresponds to the choice adopted 
in the well-known NLO analysis of the parton distribution functions by Gliick, Reya and 
Vogt. The difference between the two forms of NLO evolution solutions is fairly small, in 
this case. One might suspect that not only the scale Q 2 ni = 0.16 GeV 2 but also the scale 
Q\ n . = 0.34 GeV 2 is not high enough for the perturbative QCD framework to be justified 
perfectly. This cannot be denied completely. Still, it is clear from our simple analysis 
that there is a qualitative difference between the two choices of the starting energy, i.e. 
Qtni = 0.16 GeV 2 and Q 2 ni = 0.34 GeV 2 . As already pointed out, the authors of [37| use an 
approximate solution of the NLO evolution equation with the choice Q 2 ni = 0.16 GeV 2 to 
estimate the tensor charges at the scale Q 2 = 0.40 GeV 2 . The reduction of the magnitude 
of tensor charge after this scale change is significant. It is about 0.75 if one uses Eq.fJMj), 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 

exact RG solution 




sv 


approximate RG solution 




1 





while it is about 0.67 if one uses Eg. (IB"?]) . (See the open squares and the crosses in FigJU) 
Undoubtedly, this enormous reduction has nothing to do with the nature of their effective 
model. It is simply a consequence of starting the NLO evolution equation at such a low 
energy scale. 

Generally, for any effective models of baryons, it is very hard to say exactly what energy 
scale the predictions of those correspond to. Probably, the best we can do at the moment 
is to follow the spirit of PDF fit by Gliick et al. [67], and use the predictions of those 



mi 



models as initial-scale distributions given at the energy scale around 600 MeV, or Qf ni 
(0.3 — 0.4) GeV 2 . In fact, such approach with use of the predictions of the CQSM has 
achieved remarkable phenomenological success for both of the unpolarized and longitudinally 



polarized PDFs 



32] 



60] . In the following, we shall therefore use the exact solution (1611) of 
the NLO evolution equation with the starting energy Q 2 ni = 0.34 GeV 2 to estimate the 
tensor charges at a desired scale from the predictions of low energy models. 

In FiglHl we compare the 1st empirical information on the tensor charges for the u- and 
(i-quarks at the renormalization scale Q 2 = 0.8 GeV 2 obtained by Anselmino et al. with the 
predictions of some low energy models as well as that of the lattice QCD. They all correspond 
to the scale Q 2 = 0.8 GeV 2 . For all the low energy models except for the covariant quark- 
diquark model of 37J, the starting energy of the evolution was taken to be Qf ni = 0.34 GeV 2 



following the discussion above. In the case of the covariant quark-diquark model, we have 
tried two choices of the starting energy, i.e. Q 2 ni = 0.16 GeV 2 and Qf ni = 0.34 GeV 2 . On 



the other hand, the predictions of the lattice QCD are given in [39(] as 



A T u = 0.857 ±0.013, A T d = - 0.212 ± 0.005, 

or 

A T g (/ ~ 1} = 1.069 ±0.018, A T g (/=0) = 0.645 ± 0.018. 

Since these predictions correspond to the renormalization scale Q 2 = 4 GeV 2 , we evolve 
those down by using eq. (1641) to obtain the corresponding values at Q 2 = 0.8 GeV 2 . One sees 
that all the theoretical predictions for the d-quark tensor charge are not largely different and 
lie within the allowed range of phenomenological extraction. On the other hand, almost all 
the theoretical predictions for A^u are larger in magnitude than the empirical one, thereby 
running off the allowed range of the empirical extraction. The prediction of the covariant 
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A T u (upper) and A T d (lower) at Q 2 = 0.80 GeV 2 
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FIG. 6: Comparison of empirical and theoretical tensor charges for the u- and d-quarks. The 1st 
column and the shaded band stand for the recent empirical determination of the u- and <i-quark 
tensor charges by Anselmino et al., corresponding to the renormalization scale Q 2 = 0.8 GeV 2 . The 
theoretical predictions shown in other columns are all transformed to the same renormalization scale 
with use of the NLO evolution equation (|64p . See the text, for more detail. 

quark-diquark model with use of the starting energy Q 2 ni = 0.16 GeV 2 is an exception. 
However, we have already pointed out a serious problem of using such a low starting energy. 

At any rate, since the choice of the starting energy for low energy models is rather 
arbitrary, one must be very careful when making a comparison between model predictions 
for the tensor charges (or more generally transversity distribution) with phenomenologically 
extracted ones. (This should be contrasted with the case of axial charges. As is widely 
known, the isovector axial charge is known to be scale independent as a consequence of 
current conservation. The isoscalar or flavor-singlet axial charge is generally scale dependent, 
for example, in the standard MS factorization scheme, because of the U^(l) anomaly of QCD 
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[7J|-[76|. However, this scale dependence is known to be fairly weak except very low energy.) 
Fortunately, we can avoid this troublesome problem of initial scale choice. The key point is 
that, since the gluon does not couple to the chiral-odd transversities, the evolutions of tensor 
charges are flavor independent. This in turn means that the ratio of two tensor charges as 
AtoI/Atu or Axq^ I=0 ^ / Axq^ 1 ^ is totally scale independent. 
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FIG. 7: Comparison of empirical and theoretical tensor charge ratio A^d/A^u, which is scale 
independent. 

Shown in FigJ7]is the empirically extracted tensor-charge ratio \Axd/ A?u\ by Anselmino 
et al. in comparison with several theoretical predictions, i.e. those of lattice QCD j^], non 
relativistic quark model (NRQM) or the MIT bag model, covariant quark-diquark model 



33] , and the CQSM [3J]. We recall that this ratio is precisely 1/4 for both of the NRQM 



and the MIT bag model. One can convince that the predictions of all the models as well 
as that of the lattice QCD are not extremely far from this reference value, although the 
prediction of the CQSM is smallest of all. Since the empirical uncertainties for this ratio is 
still fairly large, we can say that all the theoretical predictions lie within the error-bars. 

Next, in FigJBJ a similar comparison is made for the tensor charge ratio Axq^ 1 ^ / 'At<^ /=1 ^ . 
Again, the prediction of the CQSM gives the smallest value among all the theoretical pre- 
dictions. Within the large error-bars, however, all the theoretical predictions are consistent 
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FIG. 8: Comparison of empirical and theoretical tensor charge ratio AtQ^ °^ / Axq^ 1 s which is 
scale independent. 

with the phenomenological value. We emphasize once again that the tensor charge ratio as 
\Axdf Aj-u] and A^q^ 1 ^ / Aj>q( I=1 ' are exactly scale independent so that it offers a safe and 
sound basis of comparison between theoretical predictions and the empirical extractions. 
Further effort to reduce the uncertainties of phenomenological extraction would be highly 
desirable. 

As a general trend, one observes that the predictions for the tensor charge ratio 
ArQ ,( ' /=0 ' ) /Ar9^ =1 ' 1 by all the low energy models as well as by the lattice QCD are not 
extremely far from the reference value of the SU(6) quark model, i.e. 3/5. This feature of 
the tensor charges should be contrasted with that of axial charges. In FiglHl we compare the 
empirically known axial-charge ratio Aq( I=0 ^ / Aq( I=1 ^ with the predictions of several models 



and with that of lattice QCD. The empirical value here is taken from the HERM 



jS analysis 



of the longitudinally polarized structure functions of the deuteron and proton (71|. (See also 
a similar analysis by COMPASS group j3|,j2^].) One sees that fairly small empirical ratio, 
which is connected with the famous "nucleon spin crisis" , is reproduced only by the CQSM 



and the lattice QCD, while the predictions of other low energy models are more or less close 
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to that of the SU(6) quark model, i.e. A^^/A* 7 ^ = 3/5, thereby largely overestimating 
this ratio. 
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FIG. 9: Comparison of empirical and theoretical axial charge ratio Aq^ 1 °) / Aq( T ^ , which is 
approximately scale independent. 

In any case, as we have repeatedly emphasized, the possible difference between the axial 
and tensor charges of the nucleon (or more generally, the difference between the longitu- 
dinally polarized distribution functions and the transversity distribution functions of the 
nucleon) offers one of the key information for disentangling the internal spin structure of the 
nucleon. Particularly useful here, we think, is the comparison between the two ratios, i.e. 
&Tcf I=Q> ' / &T<r I=1 ' and Aq( I=0 ^ / Aq( I=1 \ As we have emphasized, the former ratio is exactly 
scale independent, while the latter has only a weak scale dependence, so that it offers a safe 
and sound basis of comparison between theoretical predictions and empirical extractions for 
them. Further effort to reduce the uncertainties of the phenomenological extraction of the 
tensor charges is highly desirable to get more definite information on the possible difference 
of these two fundamental quantities. Do we expect a spin crisis also for the transverse spins, 
or the tensor charges ? 
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IV. CONCLUSION 



In summary, we have investigated the forward limit of a particular combination of 
chiral-odd generalized parton distributions, i.e. Gt(x, 0, 0) = lim^o,*-^ [H T (x,£,t) + 
2 Ht{x, £, t) + Et{x, £, t) } as well as their lower moments within the framework of the chi- 
ral quark soliton model, with particular emphasis upon the transversity decomposition of 
the nucleon angular momentum proposed by Burkardt. We found rather strong chiral en- 
hancement near x ~ for both of the isoscalar and isovector GPDs, which reminds us of a 
similar chiral enhancement observed in the CQSM predictions for more familiar unpolarized 
distribution functions of isoscalar and isovector types. We have shown that the Gt is an 
isoscalar-dominant quantity, while the isovector component also arises as an 1/N C correction. 
In particular, from the 1st moment sum rule of Gt and Ht, we have confirmed a isoscalar 
dominance of the "anomalous tensor moments" k,t, which indicates that the Boer-Mulders 
functions for the u- and ci-quarks would have roughly equal magnitude with the same sign. 
It should be contrasted with the probable isovector dominance of the Sivers functions, or 
of the anomalous magnetic moment of the nucleon. It is therefore a very important exper- 
imental challenge to determine the relative sign and the magnitudes of the u- and <i-quark 
Boer-Mulders functions. 

We have also discussed a delicate problem, which may arise when we try to compare 
the phenomenologically extracted tensor charges with corresponding theoretical predictions. 
We emphasize that the tensor charges are strongly scale dependent quantities but the ratios 
as Axd/Axu and Arg( 7=0 ) / Ayg( /=1 ) are exactly scale independent, so that these ratios are 
expected to provide us with a safe and convenient basis of comparison between empirically 
determined tensor charges of the nucleon and corresponding theoretical predictions. 

APPENDIX A: ON THE 1ST MOMENT OF E T (x,0,0) 

In this appendix, we shall explicitly verify that the 1st moment of Et{x, 0,0) vanishes. 
At the 0(Q°) level, only the isovector part of Et(x, 0,0) survives as given by eq.( 1Ml) . Its 
1st moment can easily be written down as 

/ ~nTT~ Et =1 \x, 0, 0) dx = - — ^ H r 3 « xi (7 1 75 - * 7 2 ) \ n ) 
J-i 2M N 3 ngocc 
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= -TZ E (n\T-(xxj)\n). (Al) 

l ° n£occ 

Very interestingly, this expression resembles that of the 0(f2°) contribution to the isovector 
magnetic moment of the nucleon in the CQSM given as 

fl=i\tf) = -M^ ]T (n\r-(xx^ 7 )\n), (A2) 

raGocc 

where use has been made of the generalized spherical symmetry of the hedgehog config- 
uration. Incidentally, we already know that the time reversal invariance enforces the 1st 
moment of Et to vanish. As a consistency check of our theoretical framework, we shall ver- 
ify it explicitly in the following. The formal proof in the CQSM utilizes the invariance under 
the G5 transformation, which is a simultaneous operations of the standard time reversal and 
a flavor SU(2) rotation. In a standard representation, it is given as 

G 5 = 7 1 7 2 t 2 , (A3) 

and satisfies the following identities 

G 5 Y G^ 1 = (rf , G 5 r a G5 1 = - (r a f , (A4) 
GA{x) = G 5 HG^ = H T . (A5) 

Using these properties, it is easy to verify the relations 

(n I r • (x x 7) I n) = + {n \ t ■ (x x 7) | n), 
(n I r ■ (x x 7 7) I n) = — (n \ r ■ (x x 7 7) | n), (A6) 

where use has been made of the reality of the relevant matrix elements. These relations 
then dictates that the 1st moment of E^~^ must vanish identically, while // /=1 ) need not, 
as expected. At the O^ 1 ) level, both the isoscalar and the isovector parts of Et(x, 0,0) 
survive at the first glance. In fact, their contribution to the 1st moments take the following 
forms : 

T7 (ir) E f Z f (n\(xxj) 3 \m) (m\r 3 \n), (A7) 

—^—E ( r =1 \x,0,0)dx 
-1 2M N T v ' ' J 

N c 1 

i £3ac 77 E F - F ^ ' Ta Xl 72 I m ) ( m I T c I (A8) 

13 m£nocc,n£occ 
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However, it is not so difficult to prove that both of the above expressions vanishes owing to 
the symmetry under the G5 transformation. Since the SU(2) isospin symmetry is naturally 
respected in our effective theory this just reconfirms the general statement that the 1st 
moments Ej.~ 0,1 ^ vanish by time reversal invariance. 



28 



[1] D. Mueller, D. Robaschik, B. Geyer, F.M. Dittes, and J. Horejsi, Fortsch. Phys. 42, 101 
(1994). 

[2] F.M. Dittes, D. Mueller, D. Robaschik, B. Geyer, and J. Horejsi, Phys. Lett. B209, 325 
(1988). 

[3] X. Ji. J. Phys. G 24, 1181 (1998). 

[4] K. Goeke, M.V. Polyakov, and M. Vanderhaeghen, Prog. Part. Nucl. Phys. 47, 401 (2001). 
[5] M. Diehl, Phys. Rep. 388, 41 (2003). 

[6] A.V. Belitsky and A.V. Radyushkin, Phys. Rep. 418, 1 (2005). 
[7] P. Hoodbhoy and X. Ji, Phys. Rev. D 58, 054006 (1998). 
[8] M. Diehl, Eur. Phys. J. C19, 485 (2001). 

[9] F. Ellinghaus, W.-D. Nowak, and A.V. Vinnikov, Eur. Phys. J. C46, 729 (2006). 
[10] F. Ellinghaus. larXiv:0710T5768l (2007). 



[11] Z. Ye, [hep^x /0606061 (2006). 

[12] JLab Hall A Collaboration, M. Mazous et al., Phys. Rev. Lett. 99, 242501 (2007). 

[13] X. Ji, Phys. Rev. Lett. 78, 610 (1997). 

[14] P. Hoodbhoy, X. Ji, and W. Lu, Phys. Rev. D 59, 014013 (1999). 

[15] X. Ji, J. Tang, and P. Hoodbhoy, Phys. Rev. Lett. 76, 740 (1996). 

[16] D.Yu Ivanov, B. Pire, L. Szymanowski, and O.V. Teryaev, Phys. Lett. B550, 65 (2002). 

[17] D.Yu Ivanov, B. Pire, L. Szymanowski, and O.V. Teryaev, Phys. Part. Nucl. 35, S67 (2004). 

[18] M. Anselmino, M. Boglione, U. D'Alesio, A. Kotzinian, F. Murgia, A. Prokudin, and C. Turk, 

Phys. Rev. D 75, 054032 (2007). 

[19] A. Anselmino, M. Boglione, U. D'Alesio, E. Leader, S. Melis, and F. Murgia, arXiv: 0809 .37431 

[hep-ph] (2008). 

[20] M. Burkardt, Phys. Rev. D 72, 094020 (2005). 

[21] M. Burkardt, Phys. Lett. B639, 462 (2006). 

[22] D. Boer and P.J. Mulders, Phys. Rev. D 57, 5780 (1998). 

[23] D.W. Sivers, Phys. Rev. D 43, 261 (1991). 

[24] M. Burkardt, Nucl. Phys. A735, 185 (2004) ; Phys. Rev. D 69, 074032 (2004). 

[25] M. Burkardt, Phys. Rev. D 69, 057501 (2004). 



29 



[26] J.M. Olness, Phys. Rev. D 47 (1993) 2136 

[27] H. He and X, Ji, Phys. Rev. D bf 52 (1995) 2960. 

[28] H.-C. Kim, M.V. Polyakov, and K. Goeke, Phys. Lett. B387 (1996) 577. 

[29] H. He and X. Ji, Phys. Rev. D 54 (1996) 6897. 

[30] I. Schmidt and J. Soffer, Phys. Lett. B407 (1997) 331. 

[31] L. Gamberg, H. Reinhardt, and H. Weigel, Phys. Rev. D 58 (1998) 054014. 

[32] M. Wakamatsu and T. Kubota, Phys. Rev. D 60 (1999) 034020 

[33] P. Schweitzer, D. Urbano, M.V. Polyakov, C. Weiss, P.V. Pobylitsa, and K. Goeke, Phys. Rev. 

D 64 034013 (2001). 
[34] M. Wakamatsu, Phys. Lett. B509, 59 (2001). 

[35] A.V. Efremov, O.V. Teryaev, and P. Zavada, Phys. Rev. D 70 (2004) 054018. 
[36] M. Wakamatsu, Phys. Lett. B653, 398 (2007). 

[37] I.C. Cloet, W. Benz, and A.W. Thomas, Phys. Lett. B659, 214 (2008). 

[38] S. Aoki, M. Doi, T. Hatsuda, and Y. Kuramashi, Phys. Rev. D 56 (1997) 433. 

[39] M. Gockeler, Ph. Hagler, R. Horsley, D. Pleiter, P.E.L. Rakow, A. Scafer, G. Shierholz, and 

J.M. Zanotti, Phys. Lett. B627, 113 (2005). 
[40] B. Pasquini, M. Pincetti, and S. Bom, Phys. Rev. D 72 094029 (2005). 
[41] M. Pincetti, B. Pasquini, and S. Bom, Czech. J. Phys. 56, F229 (2006). 
[42] S. Bom and B. Pasquini, Riv. Nuovo Cim. 30, 387 (2007). 
[43] H. Dahiya and A.Mukherjee, Phys. Rev. D 77, 045032 (2008). 
[44] S. Scopetta, Phys. Rev. D 72, 117502 (2005). 

[45] V.Yu. Petrov, P.V. Pobylitsa, M.V. Polyakov, I. Boring, K. Goeke, and C. Weiss, Phys. Rev. 
D 57, 4325 (1998). 

[46] J. Ossmann, M.V. Polyakov, P. Schweitzer, D. Urbano, and K. Goeke. Phys. Rev. D 71, 
034011 (2005). 

[47] M. Wakamatsu and H. Tsujimoto, Phys. Rev. D 71, 074001 (2005). 
[48] M. Wakamatsu and Y. Nakakoji, Phys. Rev. D 74, 054006 (2006). 
[49] M. Wakamatsu and Y. Nakakoji, Phys. Rev. D 77, 074011 (2008) 

[50] K. Goeke, J. Grabis, J. Ossmann, M.V. Polyakov, P. Schweitzer, A. Silve, and D. Urbano, 

Phys. Rev. C 75, 055207 (2007). 
[51] K. Goeke, J. Grabis, J. Ossmann, M.V. Polyakov, P. Schweitzer, A. Silve, and D. Urbano, 

30 



Phys. Rev. D 75, 094021 (2007). 
[52] D.I. Diakonov, V.Yu. Petrov, and P.V. Pobylitsa, Nucl. Phys. B306, 809 (1988). 
[53] M. Wakamatsu and H. Yoshiki, Nucl. Phys. A524, 561 (1991). 

[54] D.I. Diakonov, V.Yu. Petrov, P.V. Pobylitsa, M.V. Polyakov, and C. Weiss, Nucl. Phys. B480, 
341 (1996). 

[55] D.I. Diakonov, V.Yu. Petrov, P.V. Pobylitsa, M.V. Polyakov, and C. Weiss, Phys. Rev. D56, 
4069 (1997). 

[56] H. Weigel, L. Gamberg, H. Reinhardt, Mod. Phys. Lett. All, 3021 (1996). 
[57] H. Weigel, L. Gamberg, and H. Reinhardt, Phys. Lett. B399, 287 (1997). 
[58] M. Wakamatsu and T. Kubota, Phys. Rev. D 57. 5755 (1998). 

[59] P.V. Pobylitsa, M.V. Polyakov, K. Goeke, T. Watabe, and C. Weiss, Phys. Rev. D 59, 034024 
(1999). 

[60] M. Wakamatsu, Phys. Rev. D 67, 034005 (2003). 
[61] M. Wakamatsu, Phys. Rev. D 67, 034006 (2003). 
[62] M. Wakamatsu and T. Watabe, Phys. Lett. B312, 184 (1993). 

[63] Chr.V. Christov, A. Blotz, K. Goeke, P. Pobylitsa, V.Yu. Petrov, M. Wakamatsu, and 

T. Watabe, Phys. Lett. B325, 467 (1994). 
[64] M. Wakamatsu, Prog. Theor. Phys. 95, 143 (1996). 

[65] NMC Collaboration, P. Amaudruz et al., Phys. Rev. Lett. 66, 2712 (1991). 

[66] P.V. Pobylitsa, in "Large N c QCD 2004", eds. J. Goity et al. (World Scientific, 2005), p.302. 

[67] M. Gliick, E. Reya, and A. Vogt, Z. Phys. C67, 433 (1995). 

[68] S. Kumano and M. Miyama, Phys. Rev. D 56, R2504 (1997). 

[69] A. Hayashigaki, Y. Kanazawa, and Y. Koike, Phys. Rev. D 56, 7350 (1997). 

[70] W. Vogelsang, Phys. Rev. D 57, 1886 (1998). 

[71] HERMES Collaboration, A. Airapetian et al., Phys. Rev. D 75, 012007 (2007). 

[72] COMPASS Collaboration, E.S. Ageev et al., Phys. Lett. B612, 154 (2005). 

[73] COMPASS Collaboration, V.Yu. Alexakhin et al, Phys. Lett. B647, 8 (2007). 

[74] G. Altarelli and c.G. Ross, Phys. Lett. B212, 391 (1988). 

[75] D. Carlitz, J.C. Collins, and H.A. Muller, Phys. Lett. B214, 229 (1988). 

[76] A.V. Efremov and O.V. Teryaev, JINR Research Report No. E2-88, 287 (1988). 



31 



